Twirling Elastica: Kinks, Viscous Drag, and Torsional Stress 
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Biological filaments such as DNA or bacterial flagella are typically curved in their natural states. 
To elucidate the interplay of viscous drag, twisting, and bending in the overdamped dynamics of 
such filaments, we compute the steady-state torsional stress and shape of a rotating rod with a kink. 
Drag deforms the rod, ultimately extending or folding it depending on the kink angle. For certain 
kink angles and kink locations, both states are possible at high rotation rates. The agreement 
between our macroscopic experiments and the theory is good, with no adjustable parameters. 



The coupling of viscous stresses from fluid flow to de- 
formations of elastic fibers is important in many situa- 
tions, such paper manufacture, flexible microstructures 
in MEMS devices, and the dynamics of flexible biolog- 
ical filaments. The last example includes spinning Ali- 
ments, such as bacterial flagella, DNA, and supercoiling 
colonies of mutant strains of Bacillus subtilis. The ro- 
tary motors of E. coli rotate the flagellar bundle up to 
9000 rpm Q. In DNA transcription, if rotational motion 
of the RNA polymerase is blocked, then the DNA twirls 
at typical rates of 300-600 rpm as it is pulled through 
the enzyme Finally, the mutant B. subtilis cells fail 
to separate as they divide, and form long fibers which 
rotate as they grow || . 

Although these filaments are often modeled as intrin- 
sically straight elastic rods, natural bends can have sig- 
nificant effects. For example, natural curvature in bac- 
terial flagella is crucial for generating thrust; mutants 
with straight flagella can't swim J3J. Proteins can bind 
to DNA and impose sharp, large-angle bends ||, which 
can greatly enhance the torsional stress due to viscous 
drag during transcription. Similar torsional stresses from 
intrinsic bends have been estimated to be large enough 
to affect gene activity or DNA structure even in the ab- 
sence of external anchoring || . As we shall recall below, 
kinks generically trap torsional stress in specific regions 
of a rotating filament. Finally, there is evidence that in- 
trinsic bends may arise at the hairpin loops during the 
formation of supercoiled colonies of B. subtilis 0. 

To gain intuition on how sharp intrinsic bends affect 
shape and twist in the overdamped (inertialess) regime 
of cellular motions, we study a bent elastic rod rotating 
in an extremely viscous fluid, and ask, what is the shape 
and torsional stress as a function of twirling rate? For 
simplicity we disregard Brownian effects. We first describe 
how rotation affects the shape, then we formulate and 
solve the problem using slender-body theory, and finally 
we compare our predictions with experimental results. 

The shape and stresses of a rotating elastic rod in a 
viscous fluid depend on its stress-free state. At low rota- 
tion rates, a naturally straight rod twirled about its long 
axis (which is along the z-axis, say) remains straight, but 
twists and spins about z. At higher rates this state is un- 
stable, the centerline writhes and slowly rotates about z, 



and each element of the rod rapidly spins about the lo- 
cal tangent. This motion is a hybrid of crankshaft and 
speedometer-cable motion M . 




FIG. 1. Time sequence, from left to right, showing the 
steady-state shape of a naturally L-shaped rod twirling in 
glycerol. A motor (above, not shown) rotates the rod at 200 
rpm. Gravity is along z, and x points into the page. The 
container is a Plexiglas box 31 cm x31 cm wide and 28 cm 
tall. Inset is view from below, in which the rod rotation is 
counter-clockwise. 



A naturally bent rod rotating in a viscous fluid behaves 
very differently. Consider a rod made of two straight legs 
joined at a right angle. Align one leg along the vertical 
(z-) axis and twirl it with velocity loz. Unlike the natu- 
rally straight rod, the centerline of this rod will distort 
from its unstressed state for any rotation rate, since the 
free leg experiences translational drag. This translational 
drag will wrap the free leg around z (insets, Fig. 
and twist the held leg. In comparison to the naturally 
straight rod, the torsional stress in the held leg is very 
large, for a given u>. 
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FIG. 2. Bistability at sufficiently high uj. The same rod as 
in Fig. 1, spinning at 200 rpm, but a different state, obtained 
by bending the rod upwards into a V-shape at the onset of 
twirling. Inset: view from below. 

Viscous stresses tend to straighten out a rod with a 
right-angle kink (Fig. [l|).To see why, it is convenient to 
work in the frame of the twirling rod, with an unper- 
turbed fluid velocity v = luz x r, where r=r(s, i) is the 
position of the rod at time t, and s is arclength. Suppose 
the free leg aligns parallel to the x-axis for lu = 0; i.e., 
the "L" lies in the z-x plane. For small lu ( "small" will be 
defined below), translational drag deflects the free leg in 
the y-direction, which causes the held leg to bend away 
from the rotation axis in the y-direction. Thus the held 
leg lies in a region where the flow is in the — x-direction; 
this flow deflects the held leg, rotating the joint to extend 
the rod. For small lu, the sense of the joint rotation is in- 
dependent of the exterior kink angle, a. As lu increases, 
rods with small to moderate a (L-shaped) will extend, 
whereas rods with large a (V-shaped) will fold. 

These arguments suggest there should be a critical an- 
gle dividing the ultimate (large lu) extending and folding 
behaviors, but further reflection reveals that both states 
are possible at large rotation rates. Imagine that a rod 
with a right-angle kink is bent by an external force that 
brings the two legs together so the rod looks like a "U" 
with a sharp kink. Sufficient rotational flow will cause 
folding for the same reason the V-shaped rod eventually 
folds. At high rotation rates, the viscous stresses will be 
large enough to hold the rod in this folded state even 
when the external force is released. Thus, we expect 
bistability at large rotation rates: Figs, [l] and |2| show the 
same rod and same rotation rate with different initial 
conditions. These steady-state shapes are stable against 
small-amplitude perturbations, but sufficiently large per- 
turbations lead to a transition from extended to folded 
and vice-versa. 

We now turn to the mathematical analysis. For a slen- 



der rod (radius a much less than length L) , deformations 
can be described by the configuration of an orthonormal 
material frame {ei, §2, §3} embedded in the rod |^], with 
§3 = dr/ds = r s . The rates of rotation of this frame with 
respect to arclength and time define the angular strain 
f2 and velocity uj: 

(ei) s = ft x §i, (ei){ = u)Xei, (1) 

with i = 1-3. Associated with these strains are inter- 
nal elastic stresses, which give rise to a force F(s) and a 
moment M(s) acting on the cross section at s. For an 
isotropic linearly elastic rod, the constitutive relation is 
M = Ar s x r ss + Cr s r s ■ £1, where A is the bending stiff- 
ness and C is the twisting stiffness [[|. Note M = Aft 
when r = C/A — 1. In the (overdamped) limit of zero 
Reynolds number Re, the elastic force and moment per 
unit length balance exactly with the viscous drag force 
and moment per unit length: 

F s f = 0, (2) 
M s + r s x F m = 0. (3) 

From slender-body hydrodynamics Jl(J , the leading-order 
drag force and moment per unit length are 

f = C||r s r s ■ r t + (±(r t - r s r s ■ r t ), (4) 
m = ( r r s r s • uj, (5) 

where (j- ~ 2£n ~ 47r?7/[log(L/2a) + c] (n is viscosity, c 
is a constant of order unity) and ( r w Anna 2 |l(|. Al- 
though the slender-body approximation is invalid near 
the kink where the curvature diverges, it does have the 
correct (linear) w-dependence for the translational drag 
per unit length, and will therefore give the correct (non- 
linear) scaling of the shape and torsional stress with lu. 

As emphasized in ref. kinks block speedometer- 
cable motion, since sharp curvature implies a large elastic 
energy cost for rotation about the local tangent. Thus we 
consider pure crankshaft motion of the centerline about 
the axis of rotation. Furthermore, the torsional stress 
Cr s ■ f2 due to rotational drag of the held leg is (D(( r LuL), 
much smaller than the torsional stress due to transla- 
tional drag of the free leg, which for small u> is 0(£±loL 3 ). 
As we shall verify below, the twist due to translation 
dominates that due to rotation even at large to, as long 
as L/a is sufficiently large. Hence we set £ r = for fur- 
ther simplicity. Finally, we take the kink angle a to be 
clamped at a value independent of lu. 

We desire the steady-state solution to (|l|-||) for r t — 
uj x r, with ijj — luz; as our experiments confirm there is 
no oscillatory behavior as expected at low Reynolds num- 
ber. Fifteen boundary conditions at the ends are required 
for the fifteen unknowns: r, r s , ei (§2 = r^. x ei), M, and 
F. Taking the x- and y-axes to rotate about z with rate 
lu, we demand r(0) = 0, r s (0) = z, and ei(0) = x at the 
held end. The free end experiences no force and moment: 
F(L) = and M(L) = 0. Finally, fifteen conditions 
must be enforced at the position of the joint, s = L\. 
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These conditions are the continuity of position, force, and 
moment, and r s (L^) = r s (Lj~) cosa + ei(L]~) sina, and 
ei(L^) = — r s (L^) sina + ei(L^) cos a. The twist is dis- 
continuous across the kink, since the twisting moment 
on one side can balance with a bending moment on the 
other; kinks trap torsional stress. 

We measure length in units of L, time in units of 
the bending relaxation time £j_L /A, force in units of 
A/L 2 , and moment in units of A/L. The shape of the 
rod is then controlled by the dimensionless rotation rate 
X = C±^L 4 /A (cf. |nj) and the geometrical quantities a 
and Li/L. In our macroscopic experiments of Figs. 
the rod is a steel compression spring wrapped in Teflon 
tape, immersed in glycerol, with 77 re 20.0 erg-sec/cm 3 , 
L = 29 cm, a = 0.16 cm, a re 87°, and A re C re 2.2 x 10 5 
dyne-cm 2 . For our motor speeds of 5-500 rpm, x ranges 
from 55-5500. 

For small \ the shape and stress can be obtained 
to linear order in x- For the linear calculation only, 
it is convenient to write ([l]-||) in terms of r and set 
(,± = 0i; in fact, 0i will not enter the linear calcu- 
lation. The constitutive relation, the moment balance 
equation (||) (which now has m = 0), and r s • r s — 1 
imply F = — r sss + T£lr s x r ss + Ar s , where ft = r s ■ fi, 
and A is the unknown part of r s ■ F. Force balance is thus 



Xz x r = -r, 



TQr s x r sss + (Ar s 



(6) 



Since £ r = 0, the tangential moment balance (|^) implies 
O s = 0, i.e. constant twist in each leg. 

For simplicity, consider a right- angle kink, a = n/2. If 
X = 0, then fl = A = 0, and the rod is undeformed. To 
first order in x, the deformation of both legs of the rod 
is evidently along the y-direction; therefore, 

r M ~ I sz + y{s)y < s < L x 

rW ~\Li8+(«-Li)* + y(s)? Li<a<l- {<) 

Force balance implies that A s = in both legs, y SS ss = 
for < s < L%, and x( s ~ Lf) = —y SSS s for L% < s < 1. 
The boundary conditions at the ends of the rod are 
»(0) = y,(0) = 0, and j,„(l) = y sss (l) - 0(1) = A(l) = 
0. Therefore, A = fl = for L\ < s < 1; note that the 
twist for < s < L\ is not zero: the kink has trapped 
the twist. Continuity requires y(L^) = y(Lf); also, 



x + TSlLiy el 
y — TfiLix and 

z + Vs9 &i 



z - Vs y 

y-^x + yjz, (8) 



x 



yfy 



where 



the condition on the 



ei(Lf), etc. Using 
jump in the material frame becomes TflLi = yj, con- 
tinuity of force becomes y~ ss — yf ss and A - = 0, and 
continuity of moment becomes y~ s = and Tf2 = yj s . 
Thus we find a (dimensional) twist stress of Cil — 
—(±uj(L — Li) 3 /3 in the held leg. The (dimensional) 
shape of the rod is given by 



y(s) = -^s 2 (3L 1 -s)(L 1 -L) 2 , 



and 
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s > Li. (10) 



To solve for the shape and stress when x is not small, 
we must use numerical methods. However, we can obtain 
the ^-dependence of the various quantities for )(>1 us- 
ing a simple argument. Since the rod is mostly aligned 
along the z-axis at high twirling rates (in both the fold- 
ing and extending cases), r = zs + r± with |rj_| — > as 
X 0. Ignoring the terms involving A, and all numerical 
prefactors, (^) reduces to 



1 

Yj_ssss 

X 



— z x rj 

X 
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(11) 



s < Li 



(9) 



The curvature and |rj_ | are very small except in a bound- 
ary layer near the kink at s = L%; rescaling s — L\ = 
(s — LjJ/x 1 / 4 , we find that the translational drag, bend- 
ing, and twisting forces per length all balance for large 
X if O oc x 1 ^ 4 - Cm the other hand, we can balance the 
moment CD, against the approximate translational drag 
J ds(C±w|r^|)|r^| to find |rj_| oc x~ X ^- This change in 
shape with x implies that the effective rotational friction 
coefficient Cr,eff = CQ/u) decreases as x~ 3 ^ 4 f° r large x ■ 

Equations (|l|^) are in standard form for the relaxation 
method J12j . The interval < s < 1 is replaced by a fine 
mesh, with two mesh points corresponding to s = l/2 ± . 
There is no need to introduce the intermediate tension 
variable A, which enforces fixed length, since (§) and the 
boundary condition on r s at s = imply r s ■ r s = 1 . For 
simplicity, we choose T = 1 and £\\ = Ci/2- 

Fig. H shows the total extension and torsional stress 
as a function of twirling rate for a rod with a = 87° and 
Li/L = 2/3 (as in Figs.[l]-|). As X increases, the initially 
L-shaped rod distorts and extends. At sufficiently large 
X, a new branch of solutions appears which corresponds 
to folding. The situation is reminiscent of an imperfect 
bifurcation. To make the comparison with experiment, 
we have included a net downward force per unit length 
due to gravity (the filament's linear density is 0.167 g/cm 
and glycerol has a density around p = 1.3 g/cm 3 ). There 
is no fitting in Fig. ||. The agreement between theory and 
experiment is excellent. We expect the slender-body hy- 
drodynamics approximations to work well since the cur- 
vature is gentle away from the small region near the kink 
(figs. and the two legs of the rod are separated by 

many rod radii even in the folded state (inset to fig. |^). 
The Reynolds number Re re wr max ap/r] is about 0.2 for 
co around 500 rpm (r max is the maximum transverse dis- 
placement of the rod, wr max sets the velocity scale, and 
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the rod radius a sets the length scale of the disturbance 
flow). 
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FIG. 3. (a) Theoretical predictions and experimental re- 
sults for relative ^-displacement of rod free end (s = L) versus 
dimensionless rotation rate X- (b) Theoretical prediction for 
dimensionless twist in the held leg versus \. Note the agree- 
ment with the scaling arguments for fl. a ~ 87°; Li/L = 2/3. 



now assess our neglect of £ r at large x : demand- 
ing (( r wL)/(f ds(±uj\r±\ 2 ) <C 1 amounts to requiring 
X 3/4 (a/L) 2 < 1, which is easily fulfilled for L/a « 10 3 . 
Note that (within the slender-body hydrodynamics ap- 
proximation) the twist has the x 1 ^ 4 dependence for both 
the folding and extending branches. 

Fig. H displays the "phase-diagram" for a twirling bent 
rod. Rods with large kink angles (V-shaped) fold up as 
X increases from zero, whereas rods with small kink an- 
gles extend. For intermediate kink angles, both states 
are possible at high x; there is no critical angle dividing 
extending and folding. This behavior is generic as long 
as L 1 /L is not too close to or I. 

In conclusion, we have shown how kinks affect the de- 
formation and torsional stress of a twirling elastic rod 
in a viscous fluid. Despite the linearity of the elastic 
^constitutive relations (which derive from Hooke's law) 
and the equations for viscous flow (Stokes equations), we 
find nonlinear dependences on rotation rate at high rates 
due to the change in shape. Bacterial flagella, DNA, and 
fibers of B. subtilis are sufficiently flexible that typical ro- 
tation rates can cause the extending and folding studied 
here, which could be revealed using micromanipulation. 
An important extension of our work which would be rel- 
evant to DNA would be to carry out Brownian dynamics 
simulations of a rotating flexible rod with many kinks, as 
envisioned in ref. Q . 

We thank P. Nelson for posing the question that led 
to this work, and R. Goldstein, T. Peacock, H. Stone, C. 
Wiggins, and C. Wolgemuth for discussions. We thank H. 
Stone for partial support through the Harvard MRSEC 
and the Army Research Office Grant DAA655-97-I-01 4. 
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FIG. 4. Phase diagram for twirling a rod with a kink at 
Li = 2/3L (no gravity). For x ^ 1500 and a near n/2, the 
rod is bistable, as indicated by the dashed lines. Here, 2 = 1 
is fully extended and z = 1/3 is completely folded. 



Fig. g also gives the dependence of the torsional stress 
on x, confirming our scaling arguments that O oc x 
for x f and oc x 1 / 4 for \ ^> 1. We can 
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